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A FURTHER INVESTIGATION OF- THE METEOROLOGICAL 

COALITIONS CONDUCIVE TO AIRCRAFT ICING 

By William Lewis, 1 Dwight B. KLine, S 
and Charles P. Steinmetz 


SUMMARY 

Meteorological data from flight observations in icing condi- 
tions during the winter of 1946-4-7 are presented. Data on liquid 
water content, temperature, and mean— effective drop diameter are 
shown to "be consistent with values previously proposed for maximum 
icing conditions. Data on drop— size distribution as obtained by 
the rotating— cylinder method, although consistent with measurements 
previously made, were found to be inconsistent with data on drop- 
size 'distribution derived from the ratio of the maximum diameter 
to the mean— effective diameter when the maximum diameter was 
calculated from the area of impingement on a stationary cylinder. 
The relation between temperature and maximum liquid water content 
in layer clouds is discussed and estimates are given for the 
highest values of water content to be expected in layer clouds at 
various temperatures. 


INTRODUCTION 

Over a period of several years, the NACA has conducted research 
on the prevention of ice formations on aircraft through the use of 
heat. The present phase of this research is intended to provide a 
fundamental "understanding of the process of thermal ice prevention and 
thereby promote improvement in the design of thermal ice-prevention 


1 Mr.’ Lewis, U.S. Weather Bureau, has been assigned to wort in 

collaboration with the staff of the Ames Aeronautical Laboratory 
oh~the'’NACA icing research' program. ■ *-■•■-- — 

2 Mr. Kline, U.S. Weather Bureau, has' been assigned to the Flight 
Propulsion Research Laboratory, Cleveland, also in connection 
- ’With the' NACA icing- research program,.- 
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equipment. The meteorological results of the investigation, up to 
and including the 1945—46 winter season-/ have teen presented in 
reference 1. The present report, which may be regarded as a 
supplement to reference 1, presents the results of additional 
flight measurements made during the winter of 1946-47. These data 
include results of flight observations with the XB— 2i<M and XB— 25® 
airplanes .operated by the Flight Propulsion Research Laboratory at 
Cleveland, as well as the C— 46 airplane operated by the Ames 
Aeronautical Laboratory, Moffett Field, ‘ Calif . 

Appreciation for- cooperation in this investigation is extended 
to United Air Lines, Inc., the U.S. Weather Bureau, and the Air 
Materiel Command of the Army Air Forces. In particular. Major 
James Murray, Air Materiel Command, and Mr. Lyle Reynolds, United 
Air Lines, who were pilot and copilot, respectively, of the 
C— 46 airplane, contributed materially to the research program. 


APPARATUS AND METHOD 

The equipment and test methods employed in the research of this 
' report were identical to those described in reference 1 with the 
exception of a few changes discussed in the following paragraphs. 



Rotating— Cylinder Measurements 

The values of liquid water content and mean-effective drop 
diameter 3 presented herein were all calculated from the amounts 
of ice collected on four rotating cylinders, l/8, l/2, 1—1/4, and 
3 inches in diameter. The assembly used on the C— 46 airplane is 
shown in figure 1* A similar apparatus was used by the Flight 
Propulsion Research Laboratory. 

Previous calculations of liquid water content and mean- 
effective diameter frcm rotating— cylinder data obtained in flight 
■'have been based upon the true airspeed of. the airplane. . (See 
'■’references 1 and 2.) This procedure involves the tacit assumption 
that the local velocity at the point where the cylinders are 
'exposed is equal to the true, airspeed of the airplane. In order 
to check- the validity of this .assumption for the- installation 


3 Defined in the appendix of reference 1 as the volume median 

dionfeter having the property -that there, is as much water in the 
drops larger than the. volume., median diajtjeter as there is in 
drops smaller than, :th©- -volume. median diameter. * 
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on the C— 46 airplane , the local velocities at the points of exposure 
of- the cylinders were measured for a range of airspeeds. It was 
found that the local velocity was approximately 12 percent higher 
than the true airspeed over the entire region in which the cylinders 
were exposed. The values of mean-effective diameter and liquid 
water content presented in this report, calculated from observations 
on the C— 46 airplane, are based upon the local velocity. The data 
presented herein from the Flight Propulsion Research Laboratory are 
based on true airspeed, since local velocity data were not available 
■for those installations. 


Area of Drop— Impingement Measurements 

The apparatus for the measurement of the area of drop impinge- 
ment as used on the C— 46 airplane is shown in figure 2. This device 
consists of a cylinder 5 inches in diameter on which the angle from 
the stagnation point was marked in intervals of 10°. A means was 
provided for removing the ice accretion by rotating the cylinder 
against a scraper which was mounted directly behind, the cylinder. 

The angle between the aft edge of the ice formations and the 
cylinder stagnation point- was read visually to the nearest 5°. 

The ice formation was not allowed to became large enough to signifi- 
cantly modify the circular cross section of the cylinder. 


Icixg Rate Meter 

An icing rate meter of tha rotating disk type was used on the 
airplanes operated by the Flight Propulsion Research Laboratory. 
This instrument was similar in principle to the one described in 
reference 2. Data from this instrument are presented as icing rate 
in inches per hour collected on the edge of the disk. The data are 
not presented in terms -of liquid water content, since the density 
of the ice and the collection efficiency of the disk are not known 
with sufficient accuracy. 


RESULTS AND DISCUSSION 

The data from the 1946-47 observations have been prepared 
in a form similar to that used for presentation in reference 1. 
Tables I and H present a summary of the data for both laboratories 
for all flights during which icing conditions were measured. 

Figure 3 presents the relation between liquid water content and 
mean-effective drop diameter and figure 4 shows liquid water content 
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as a function of free-air temperature.' .The curves denoting intensity 
of icing in figure 3 represent the rate, of ice accretion on a 3— inch 
diameter cylinder at 200 miles per hour as specified hy the 
IT.S, Weather Bureau for reporting icing intensity from mountain-top 
observation stations. 

All of the values of liquid water content and meanr-effective 
drop diameter listed in tables I and II are -within the range of 
values observed during the previous season except the largest value 
of drop size from flight 72. In this case the collection efficiency 
was the same, within the limits of errors of measurement, for each 
of the four cylinders.' This corresponds to a value of moan- 
effectivo diameter of at least 150 microns. The liquid water 
content was 0.04 gram per. cubic meter in this case. 


• Icing Conditions in Altos tratus Clouds 

■ In contrast to the experience of the 1945--46 season in which 
nearly all of the. altos tratus clouds observed were composed of ice 
crystals, a largo altostratus cloud Bystem composed mostly of wator 
drops was 'encountered on flight lOg in the zone of convergence ahead 
of a low-pressure area. This cloud was formed in tropical marine 
air over Louisiana at a temperature only slightly below freezing. 
Altostratus clouds containing water drops were also encountered 
during four flights by the flight Propulsion Research Laboratory. 
Three of those were Just outside cyclonic precipitation areas and 
the fourth was ahead of a cold front. A more thorough investiga- 
tion - of the structure of clouds Associated with, fronts and low- 
pressure systems will be required to determine the most probable 
location and extent of icing conditions in such cloud systems. 

The following summary of the icing characteristics of alto- 
stratus clouds as compared with stratocumulus clouds in the same 
general area has been prepared from data obtained by the Plight 
Propulsion Research Laboratory. 




Cloud type 


Number of runs 


Altostratue 


average * 0.19 

Liquid water content median i .18 


range \ .12 to .30 


Mean-effective 
drop diameter 


Temperature 



Strat ocumulus 


0.21 

.18 

.06 to .50 


13 microns 
12 microns 
7 to 36 microns 


12 ° F 
15 ° F 

-11 to 28 ° F 


Although, insufficient data Are presented here to permit definite 
conclusions to be drawn, it ia noted tliat there 1 b a tendency for 
altostratue clouds to have larger drops and more uniform conditions 
than stratocumulus clouds. 


Re-examination of Previously Proposed 
Design Conditions 

Tentative estimates of the most severe icing conditions likely 
to be encountered in the course of all-weather transport operations 
in the United States were presented in reference 1 to serve as a 
guide in the design of ice— prevention equipment These estimates, 
which were based on data obtained during the 19^5— U6 season and 
before are given below: 
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. It is seen from an examination of the data in figure 3 that 
these maximum icing conditions were not equaled or exceeded during 
the 19L6-47 observations. The general range and frequency of 
values of liquid water content and mean-effective diameter in layer 
clouds are very similar to the results from the previous season. 

Only a small amount of data were taken in. cumulus clouds in 1946 - 47 . 
These observations all fall within the range established by the 
19^5-46 observations for cumulus clouds. • 


The Relation Between Maximum Liquid Water Content 
and Temperature' in- Layer Clouds 

It was pointed out in reference l.that insitfficient data were 
available from layer clouds at low temperatures to provide the 
basis for an estimate of the relation, botwo.cn temperature and 
maximum liquid water content in layer-type clouds. Th*. data 
presented herein, while still rather scanty, include observations 
from 12 flights in layer clouds at temperatures of 10° F or lower 
and two flights below —10° F. These are. believed to provide a 
sufficient basis for a tentative estimate of maximum liquid water 
content as a function of temperature in layer clouds. 

It was suggested In reference 3 that the maximum liquid water 
content likely to occur in stratus clouds Is the amount that would 
be produced by adiabatic lifting through an Interval of 3000 feet 
above the condensation lev-1. Subsequent experience, and improved 
methods of measuring liquid water content indicate that tho actual 
water content is generally substantially loss than tho theoretical 
value. The estimate of 3000 feet as the maximum thickness 'of a 
continuous etratocumulus or stratus layer appears to be approxi- 
mately correct. In figure 4 , curve A represents one-half of tho 
liquid water content which would be obtained by adiabatic lifting 
through a pressure altitude interval of 3000 feet from the condensa- 
tion level. This curve falls very close to’ the points representing 
the highest observed values of liquid water oontent in layer— type 
clouds. Since it is reasonable to expect that a larger sample of 
data would Include higher values pf water content, the curve B 
(fig. 4 ) , which represents two-thirds of the liquid water content 
produced by adiabatic lifting through 3000 feet, is proposed as an 
estimate of the highest values of liquid water content to be 
expected In layer clouds,: This, curve indicates a maximum liquid 

water content of 0.8 gram per cubic meter at 20° F which is in 
agreement with the estimate given in reference 1. The maximum 
liquid water content for lower temperatures is 0.5 gram por cubic 
meter at 0° F and. 0.25 gram per cubic meter at -20° F. 
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Typical Ioln& Conditione 

Estimates of the most severe Icing conditions likely to be 
encountered in the course of all -weather transport operations in 
the United States have been presented in the foregoing section. 

Bata on typical or average icing conditions and on the relative 
frequency of various values of liquid water content and drop size 
are also of interest. 

The highest values of liquid water content measured during each 
of 21 flights in cumulus clouds anc. 51 flights in layer— type clouds 
are presented in figure 5 in the form of ogives (cumulative frequency 
curves, reference- k) . These curves include data from the Ames 
Aeronautical Laboratory for 19^5— 46 and 1946-47 and from the Flight 
Propulsion Research Laboratory for 1946-47. The median value of 
maximum liquid water content par flight is 0.7 6 gram per cubic meter 
for cumulus clouds and 0.28 gram per cubic meter for layer clouds. 
These values approximate those given in reference 1 for typical 
icing conditions. It is also noted from figure 5 that- 90 percent 
of the flights in cumulus clouds encountered less than 1.2 grams 
per cubic meter and 90 percent of the flights in layer clouds 
encountered less than 0.5 gram per cubic meter. 

Ogives plotted from observations of mean-effective diameter ■ 
and maximum diameter made during the 1946-47 season are presented 
in figure 6. It is noted that 50 percent of the observations of 
mean-effective diameter fall in the relatively narrow range from 
11.2 to 16.2 microns, and 90 percent are Iobs than 22 microns. 

Fifty percent of the observations of maximum diameter are between 
12.6 and 20 microns and 90 percent are below 28 microns. 


These distribution curves indicate that the icing conditions 
most commonly encountered are much less severe than the estimated 
maximum conditions. Thus, if it were assumed that most cumulus 
clouds and the most severe 10 percent of icing conditions in layer • 
clouds could be avoided by proper meteorological navigation, it 
would only be necessary to protect against 0,5 gram per cubic meter 
at 13 microns or 0.3 gram per cubic metey at. 20 microns. The extent 
to which meteorological navigation can be relied upon, however, can 
only be determined by an extensive study of the distribution of 
icing conditions in various weather situations and an analysis of 
air traffic control procedures. 
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Maximum Drop Size and Drop-Size. Distribution 

In the" analysis of the test data the fact us noted that, in 
many cases, the maximum drop diameter as calculated from the area 
of impingement on the fixed cylinder was equal to or only slightly 
exceeded the corresponding .mean-effective diameter calculated from 
the rotating— cylinder data.’ In a few cases the indicated maximum 
diameter Was less than. the indicated mean— effective diameter. This 
would indicate that in a majority of cases the size distribution was 
fairly uniform. The size distribution obtained by the rotating— 
cylinder method, on the Other hand, frequently indicated broad 
distributions in cases where the comparison of the mean and the 
maximum diameters indicated uuif ohm. drop size. In order to check 
the consistency of. these two methods- of measuring drop-size 
distribution, they were’ expressed in terms of a common scale. To 
do this the assumption was made that’ the value of drop diameter 
contributing 10 percent of the water content in the assumed size 
distributions B, C, D, and E (reference &) corresponds to- the 'value 
of maximum drop diameter derived from the area of- impingement on 
the stationary cylinder. On the basis of this assumption,' the ratio 
of the maximum diameter as measured by the area of impingement method 
to the. mean— effective diameter as measured by the rotating— cylinder 
method was used to define a scale of size distributions as follows: 


Size distribution designation Maximum diameter 
( defined in reference 2) mean-effective diameter 


A 

B 

C 

D 

E 


below 1.16 
,1.17 to 1.1(1 
1.42 to 1.62 
1.63 to 1.87 
1,88 and over 


A comparison of - the drop— size distribution obtained by the two 
methods is shown in the following frequency table: 
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Hurnbor of observations 

Size distribution determined hy the ratio of 


maximum diameter to mean— effective, diameter 


Size distribution by 
the rotating-cylinder 
method 


A 

3 

C 

D 

■3 

Total 

A 

21 

10 

1 

1 

1 

34 

B 

3 

6 ’ 

3 

1 

1 ' 

14 

0 

8 

b 

0 

1 

2 

15 

x> 

4 

b 

0 

0 

0 

8 

E 

8 

b 

2 

5 

2 

21 

Total 

44 

28 

6 

8 

6 

92 


It is seen from the foregoing frequency table that the data on 
drop— size distribution determined by the ratio of maximum diameter 
to mean— effective diameter indicate a preponderance of fairly uniform 
drop— size distributions; whereas the data on drop— size distribution 
obtained from the rotating— cylinder method indicate a larger number 
of very nonuniform distributions. - 

The correlation coefficient showing the degree of agreement 
between the results of the two methods of measuring size distribution 
was computed from the foregoing table by Pearson's product moment 
formula (reference 4) and found to be 0.19. This low correlation 
between the results by the two methods indicates that one or both 
methods must be regarded as unreliable and that therefore the informa- 
tion presented herein on drop— size distribution must, at present, bo 
regarded with scsne skepticism. The values of mean-effective diameter 
presented are nevertheless regarded as being fairly accurate. 

An Inspection of the data in the frequency table shows that tho 
agreement in drop— size distribution would not be materially improved 
by a modification of the assumption used to reduce the data to a 
common scale. For example, if the scale had been chosen to give 
" E n distribution for a larger fraction of the maximum drop— size 
data, the improved agreement in the lower portion of the table 
would be offset by corresponding, changes, in the upper portion. 

One possible explanation for this discrepancy in the determina- 
tion of drop— size distribution lies in the effect of tho acceleration 
in the flow of air around the fuselage in locally modifying the water 
content, drop— size distribution, and velocity at the points where ' 
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the rotating cylinders are exposed. Since the l/8-lnch cylinder 
is exposed nearly twice as . far from the aide of the fuselage, as the 
3— inch cylinder, it would appear that, the local effects might apply 
variously to the different -cylinders, thus giving rise to a false 
relationship "between cylinder diameter and relative collection 
efficiency. Only a small change in the curvature of the line 
defining this relationship is sufficient to produce a significant 
change in the indicated drop-size distribution. This effect could 
he measured by the exposure of four rotating’ cylinders of ejual 
diameters at positions normally occupied by the rotating cylinders. 
Any differences in the amounts of ice collected would be due to tho 
local acceleration effects Just mentioned. 

Another possible explanation of the discrepancy is the 
possibility that flow around the ends of the stationary cylinder 
caused the observed width of the area of impingement to be less 1 

than would occur on a cylinder of infinite longth, This effect is 
believed to be unimportant, however, since the edges of the ice 
formations were observed to be straight and parallel to' the axis 
of the cylinder. 

The frequency of various values of maximum drop diameter are 
presented in figure 6. It is seen' from this curve that the maximum 
diameter was less than 20 microns in 75 percent of observations and 
less than 3.0 microns in 93 percent. In a comparison of the two 
curves of figure 6, it should bo remembered that the curve for 
maximum diameter is based upon a much smaller Bample of data than 
tho curve for mean-effective diameter. 


A Further Check of the Icing Intensity Scale Proposed 
in Reference 1 for Forecasting Purposes 

It was pointed out in reference 1 that, while fairly Tellable 
estimates of the liquid-water content in clouds can be made, the 
size of the drops remains unpredictable. For this reason, a scale 
of icing intensity based upon liquid water content alone was 
proposed as an aid in the preparation of icing forecasts. This 
scale was found to agree with the icing intensity scale used by 
the Weather Bureau in 7& percent of the l^b^—k6 observations. The 
I 9 I 16 — b 7 observations have been used to check the general validity 
of the proposed scale, since these data are independent of those 
used to define the scale. The following table presents the 19^6 — hj 
data in the same form used in reference 1 for the 19^5— 46 data. 
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; Cloud type 

l 

t 

1 



. Range of 
water 
content 

1 i 

Number of observed cases of icing of 1 
various intensities 1 

Alternate j 
scale of l 

Weather Bureau scale cf 
‘ intensity 

icing s 

icing [ 

intensity ; 

i 

Trace ‘ 

Light- 

Moderate 

■ Heavy 

1 

! Layer 
; clouds 
i 

0-0. ii 
■ 0.12-0.68 
0.69-1.33 
over I.33 

Trace : 
Light ■ 
-Moderate j- 
Heavy 

. U : 

0 ! 
0 

7 ' 
79 ", 
0 * 

0 - 

O OJ O O 

0000 

1 

1 

0-0.07 

Trace ! 

. 0 ! 

0 . 

0 

f *■ i 

i 0 1 

1 Cumulus 

0.08-0.49 

Light [ 

3 « 

13 

0 

; 0 

1 clouds 

• 0. $0-1.00 ; 

Moderate-- i 

0 1 

2 • 

3 

j. 0, ; 

_ 

1 over 1.00 

.. 

Heavy 

0 ; 

O 

0 

Lli 


This table shows agreement in 76 percent of the observations. 

The alternate scale indicates icing intensity one degree higher in 
18 percent of the cases and one degree lower in 6 percent. The corre- 
sponding figures for the 194 5-46 data were 78 percent, 17 percent and 
5 percent, respectively. Thus, the agreement is nearly as good for 
the independent data, which indicate, in general, that approximately 
this degree of agreement can be oxpected. 

CONCLUDING REMAKES 

In addition to verifying the specifications of maximum icing 
conditions proposed previously, the data presented herein have been 
used to define the relation between temperature and maximum liquid 
water content in layer-type clouds as follows: 

T em perature Maximum liquid water content 
in layer clouds 

20 ° F 0.8 gm/m 3 

0 ° F .5 gm/m 3 

- 20 ° F .25 gnyW 3 

Data on drop-size distribution as obtained by the rotating— cylinder 
method, although consistent with measurements previously made, were 
found to be inconsistent with data on drop-size distribution derived 
from the ratio of the maximum diameter to the mean— effective diameter 
when the maximum diameter wets calculated from the area of impingement 
on a stationary cylinder. In spite of the inconsistency it is 
believed that the data on mean-effective diameters are fairly reliable. 

Ames Aeronautical Laboratory,. 

National Advisory Committee for Aeronautics, 

Moffett Field, Calif. 
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Figure 1.- Rotating cylinder apparatus used on the C-46 airplane 
during icing research in the 1946-47 winter. 
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Apparatus used to measure the area of drop impingement 
installed on the C-46 airplane. 


Figure 2.- 




liquid i/v'ATCQ CONTENT, ^An 3 







FIGURE 4. - UQUIO WATER. CONTEKT AS RELATED TO TElAPEftfVTURC 
.IN \C|NG C.LOOOS. 
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FIGURE 5.- OGIVES OF KiAVlMUM 
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FIGURE 6 - OGIVES OF MEAN EFFECTIVE DlAMCTER AND NVtf IMUt'T DIAMETER 
OF DROPS l N IC/NG CLOUDS. 
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